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Abstract. Quantitative spectroscopic analyses of two 
faint blue stars (V=19™5 - 20™0) in the globular 
cluster M 15 are presented. Their derived T e g, gravi- 
ties and absolute magnitudes (T e g =24000K, logg=5.2, 
M v =4^3; T cff =36000K, logg=5.9, M v =4'°7, respec- 
tively) are matched very well by models for the Extreme 
Horizontal Branch (EHB). Both stars are bona- fide subd- 
warf B stars making M 15 only the second globular clus- 
ter (after NGC 6752) for which the existence of sdB stars 
has been proven spectroscopically. While the helium abun- 
dance (one tenth solar) of Fl-1 is typical for sdB stars, 
F2-2 surprisingly turned out to be a helium rich star, the 
first to be reported as a member of a globular cluster. In 
the field population of the Milky Way such stars are rare 
(less than 5% of all sdB stars). From its proximity to the 
helium main sequence, it is speculated that F2-2 may be 
a naked helium core, i.e. an Extreme Horizontal Branch 
star which lost (almost) all of its hydrogen-rich envelope. 



Horizontal Branch (EHB) stars. Such objects are now be- 
lieved to be the dominant source of UV radiation caus- 
ing the UV upturn phenomenon in elliptical galaxies and 
galaxy bulges (Dorman et al., 1995| ). While there are more 
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than a thousand sdBs known in the field of our galaxy 
(Kilkenny et al., 1988), only in one globular cluster bona- 
fidc sdB stars were show n to e xist (NGC 6752, Heber et 
ah, |l986j , Moehler et al., |l996|) . Several claims that sdB 
stars had been found in other globular clusters too could 
not be confi rmed by spectr oscop ic analyses. Moehler et al. 
(|!995i |1996| , de Boer et al. |1995| , hereafter Paper I,III, and 
II respectively) found that all "classical" BHB stars (i.e. 
stars with T e g < 20000 K and logg < 5) in several clusters 
exhibited too low masses compared to standard evolution- 
ary theory, whereas the sdB stars' masses in NGC 6752 
were in good agreement with the canonical mass of 0.5 
M . 

It is therefore of great importance to find and anal- 
yse sdB stars in other globular clusters. In this letter we 
present follow-up spectroscopy and spectral analyses of 
faint blue stars (19 m < V < 20 m and -0™28 < (B-V) < 
-0™12) in the globular cluster M 15, which were discovered 
recently by Durrell & Harris ( 1993 ). 



1. Introduction 



2. Observations and Reduction 



Subluminous B stars (sdB) form the extreme blue end 
(T Gff > 20000 K and logg> 5) of the Horizontal Branch 
(Heber et al., 1984) and are therefore termed Extreme 
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Two of the four candidates (F2-1 and F2-3) could not be 
observed reliably from the ground due to nearby red neigh- 
bours (see Table 1). The remaining candidate stars were 
observed with the focal reducer of the 3.5m telescope at 
the Calar Alto observatory using grism #3 (134 A/ mm 
and a slit width of l'.'S, resulting in medium resolution 
spectra covering the wavelength range 3250 - 6350 A. We 
also obtained low resolution spectrophotometric data, us- 
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ing a slit width of 5" and binning the spectra by a factor 
of 2 along the dispersion axis. 



Number 


V 


B-V 


Q2000 


#2000 


Fl-1 


19™468 


-0*?160 


21 h 30 m 12 i !3 


+12°03'45'.'6 


F2-1 


19^101 


-0^126 


21 h 29 m 39 ! !6 


+12°07'26'.'3 




18™555 


+1™531 


21 h 29 m 39 ! !6 


+12°07'26'.'8 


F2-2 


19™983 


-0*?231 


21 h 29 m 34=7 


+12°09'19'.'l 


F2-3 


19™956 


-0" 1 276 


21 h 29 m 46 s 6 


+12°06'53'.'7 




18™027 


+0™702 


21 h 29 m 46 i :6 


+12°06'54'.'3 



Table 1. Positions and magnitudes of the sdB candidates and 
their close neighbours 



Since there is no built-in wavelength calibration lamp 
we observed wavelength calibration spectra only at the 
beginning of the night. The observation and calibration 
of bias, dark current and flat-field was performed as de- 
scribed in paper I and III. As the focal reducer produces 
rather strong distortions we extracted that part of each 
observation that contained the desired spectrum and some 
sky and straightened it, using a program written by Dr. 
O. Stahl (priv. comm.). We applied the same procedure 
to the wavelength calibration frames. Thereby we could 
perform a good two-dimensional wavelength calibration 
(using the MIDAS Long context) and sky-subtraction (as 
described in paper I). Correction of atmospheric and in- 
terstellar extinction as well as flux calibration was again 
performed as described in paper I using the flux tables of 
Massey ( |1988| ) for BD+28°4211. 

3. Analyses 

3.1. Fl-1 

The spectra show broad Balmer lines and a weak He I 
line at 4471 A typical for sdB stars. The low resolution 
data unfortunately showed some excess flux towards the 
red, probably caused by a red star (V = 19™99, B— V = 
+0™55) 3" away. Fitting only the Balmer jump and the B 
and V photometry of Durrell & Harris (1993) we get an 
effective temperature of 24000 K (cf. paper III). We used 
the helium- and metal-poor model atmospheres of Heber 
(cf. Paper III) to analyse the medium resolution spectrum 
of Fl-1. As the resolution of the spectrum varies with 
wavelength we convolved the model atmospheres with dif- 
ferent Gaussian profiles for the three Balmer lines. The 
appropriate FWHM were determined from the calibra- 
tion lines close to the position of the Balmer lines. We 
used FWHM of 8.7 A for H 5 , 7.1 A for H 7 , and 9.8 A for 
Hg. We fitted each line separately and derived the mean 
surface gravity by calculating the weighted mean of the 
individual results. The weights were derived from the fit 
residuals of the individual lines. We thus get a log g value 
of 5.2±0.14 for a fixed effective temperature of 24000 K. 



The fits to the Balmer lines are shown in Fig. 1. Note in 
passing that higher temperature and gravity results if we 
ignore the Balmer jump and derive T e g and logg simul- 
taneously from the three Balmer line profiles alone: The 
smallest overall residual is achieved for T c s — 29700 K and 
logg = 5.9. These values, however, are inconsistent with 
the low resolution data. We already noted similar incon- 
sistencies between Balmer jump and Balmer line profiles 
in paper I, which are probably caused by insufficient S/N 
in the Balmer line profiles. The Hei 4471 A line is consis- 
tent with a helium abundance of about 0.1 solar. Using 
the routines of R. Saffer (Saffer et al., 1994 ) and a fixed 
mean FWHM of 8 A we get internal errors of 600 K and 
0.13 dex, respectively. We take the standard deviation of 
log g as external errors and assume an external error of T c ff 
of ± 2000 K (cf. Paper III) . Using the same method as in 
Papers I and III we get a logarithmic mass of — 0.423±0.20 
dex, corresponding to (0.38^0.14) M Q . For a detailed dis- 
cussion of the errors entering into the mass determination 
see Paper III. 



T3 




-50 

AA (A) 

Fig. 1. The best fitting models for T cfi = 24000 K compared 
to the Balmer lines of of Fl-1 (logg = 5.0 (H 4 ), 5.3 (H 7 , Hp)). 
The tickmarks mark steps of 10% in intensity. 



3.2. F2-2 

Unlike Fl-1 and other sdBs the spectrum of F2-2 is not 
dominated by Balmer lines but displays a prominent He 1 
absorption line spectrum in addition to the Balmer lines 
(see Fig. 2). The Balmer lines are much weaker and nar- 
rower than in Fl-1. In Fig. 2 we compare the spectrum of 
F2-2 (smoothed by a box filter of 5 A width) to t hat o f 
a He-sdB in the field (HS 0315+3355, Heber et al. ~ 



1996|) , 



a spectral type that describes rare helium-rich variants 
of the sdB stars (Moehler et al., 1990| , see also below), 
The helium line spectra of both stars are very similar, 
while the Balmer lines in HS 0315+3355 are considerably 
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weaker than in F2-2 (see Fig. 2), indicating an even lower 
hydrogen content of the former. 




4000 



4500 
A[A] 



5000 



Fig. 2. The medium resolved spectrum of F2-2 and the iden- 
tified absorption lines. To allow a better identification of the 
lines the spectrum has been smoothed with a box filter of 5A 
width. In comparison (thick line) we show the spectrum of a 
helium-rich sdB star (HS 0315+3355) in the field that has been 
convolved to the same resolution. The tickmarks mark steps of 
10% in intensity. 



The admittedly rather low S/N unfortunately allows 
only a coarse analysis. Since an analysis of individual spec- 
tral lines is impossible we choose to fit the entire useful 
portion of the spectrum (4000A - 5000A) simultaneously. 
It turned out that F2-2 is somewhat hotter than Fl-1 and 
we therefore used the updated grid of NLTE model atmo- 
spheres calculated by S. Dreizler (see Dreizler et al. 1990 ). 
The model atmospheres include detailed H and He model 
atoms and the blanketing effects of their spectral lines but 
no metals. Using Saffer's fitting program T e g, logg, and 
helium abundance were determined simultaneously from 
a small subgrid (T cff = 35, 40, 45kK; logg = 5.5, 6.0, 6.5; 
He/ (H+He) = 0.5, 0.75, 0.91, 0.99, by number). 

Since the Hei 4144A line is not included in the models, 
this line is excluded from the fit. 

An effective temperature of 36000 K, logg of 5.9, and 
a helium abundance of N/f- e /(N#+N# e ) = 0.87 (by num- 
ber) resulted. The fit to the spectrum is displayed in Fig. 
3. Since the noise level is rather high it is difficult to esti- 
mate the error ranges of the atmospheric parameters. As 
a test we changed the fitted spectral range as well as the 
placement of the continuum, which resulted in quite small 
changes of T ofr (« 2000K) and logg (« 0.2 dex). The he- 
lium abundance Nff e /(N#+Nij e ) ranged from 0.5 to 0.9 
by number, indicating that helium is overabundant by at 
least a factor 5 with respect to the sun. As conservative 
error estimates we adopted ± 4000 K and 0.5 dex for the 
errors in T c g and logg. 
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Fig. 3. The spectrum of F2-2 and the best fitting model (thick 
line, r cff = 36000 K, logg= 5.9, N He /(N H +N He ) = 0.87). To 
allow a better identification of the lines the observed spectrum 
has been smoothed with a box filter of 5A width. The tickmarks 
mark steps of 10% in intensity. 



3.3. Radial velocities 

For Fl-1 we derive a mean heliocentric velocity of 
— 142 km/sec from the Balmer lines with an estimated er- 
ror of about ±40 km/sec (due to the mediocre resolution, 
the limited accuracy of the wavelength calibration, and the 
shallow lines). This value is within the error limits consis- 
tent w ith th e cluster velocity of -107.09 km/sec (Peterson 
et al., |l989 ). For F2-2 we could not derive reliable radial 
velocities from single lines due to the low S/N of the spec- 
trum. Instead we cross correlated the normalized spectrum 
with that of the field He-sdB HS 0315+3355 (convolved to 
the same resolution). This procedure resulted in a helio- 
centric velocity of ~ —70 km/sec, which - due to the large 
errors - does not contradict a cluster membership for F2-2. 



4. Discussion 

The spectroscopic analyses of two faint blue stars in the 
globular cluster M 15 show that both stars are bona-fide 
subdwarf B stars. In Fig. 4 and 5 we compare their posi- 
tions in the (T e g, logg)- and the (T e g, My)-diagrams to 
those of EHB stars in NGC 6752 (from paper III) as well 
as to the predictions from Horizontal Branch stellar mod- 
els (from Dorman et al. 1993). As can be seen from this 



comparison, their evolutionary status is well described by 
models for the Extreme Horizontal Branch (EHB). Hence 
M 15 is only the second globular cluster (after NGC 6752) 
for which the existence of EHB stars has been proven spec- 
troscopically. While the helium abundance of Fl-1 is typ- 
ical for sdB stars (i.e. subsolar), F2-2 surprisingly turned 
out to be a helium rich star. 

This is the first time ever that a helium rich sdB star 
has been reported in a globular cluster. In the field of the 
Milky Wa y only 5% of the sdB stars are helium-rich. Jef- 
fery et al. ( 1996|) list 48 such stars, while the latest (unpub- 
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□ M 15 (Paper I) 
H NGC 6752 (Paper III) 
■ this paper 




□ M 15 (Paper 1) 
H NCC 6752 (Paper III) 
■ This paper 



Fig. 4. The physical parameters of the sdBs in M15 compared 
to the results of Papers I and III and to theoretical expec- 
tations. The solid lines are the Zero Age HB and the Termi- 
nal Age (core h elium exhaustion) HB for [Fe/H] = —2.26 of 
Dorman et al. (1992). The short dashed line gives the posi- 
tion of the helium main sequence (Paczynski 1971). The long 
dashed-short dashe d line s give post-EHB evolutionary tracks 
by Dorman et al. (1993), labeled with the total mass of the 
EHB star. 



lished ) vers ion of the catalog of hot subdwarfs (Kilkenny 
et al., 1988] ) lists more than 1,000 hydrogen-rich sdBs. 

The helium-rich sdB has an absolute visual brightness 
of about 4™7, which places it at the very faint blue end 
of the EHB as seen in the colour-magnitude diagram of 
NGC 6752. F2-2 may even be hotter than any EHB star 
in NGC 6752. From its proximity to the helium main se- 
quence in Fig. 4 and 5 it might be tempting to regard 
F2-2 as a naked helium core, i.e. as an Extreme Horizon- 
tal Branch star which lost (almost) all of its hydrogen-rich 
envelope. 

Why didn't we find any helium-rich sdBs in NGC 6752? 
All the EHB stars that have been analysed in NGC 6752 
(inclu ding the three faintest ones seen in Buonanno et al. 



1986 ) are helium-poor sdB stars (Paper III). This could 
either mean that there are no helium-rich subdwarfs in 
NGC 6752 or that t hey ar e just below the detection limit 
of Buonanno et al. (1986). One should certainly keep an 
eye on newer and deeper CMDs of globular clusters to see 
whether other He-sdB candidates show up. 
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Fig. 5. The absolute V magnitudes and effective temperatures 
as g iven above compared to theoretical tracks by Dorman et 
al. (1993, details see Fig. 4). Also shown are the data for the 
stars analysed in papers I and III. 
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